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LEAST-INDEX  RESOLUTION  OF  DEGENERACY 
IN  LINEAR  COMPLEMENTARITY  PROBLEMS 
WITH  SUFFICIENT  MATRICES 

by  Richard  W.  COTTLE  and  Yow-Yieh  CHANG 
ABSTRACT 

This  paper  deals  with  the  Principal  Pivoting  Method  (PPM)  for  the  Linear  Com¬ 
plementarity  Problem  (LCP).  It  is  showm-herA-that  when  the  matrix  M  of  the 
LCP  {q,M)  is  (row  and  column)  sufficient,  the  incorporation  of  a  least-index 
pivot  selection  rule  in  the  PPM  makes  it  a  finite  algorithm  even  when  the  LCP 
is  degenerate. 


1.  Introduction 

The  Principal  Pivoting  Method  (PPM)  for  the  Linear  Complementarity  Problem  (LCP) 


w  —  q  +  Mz 

(1) 

tv  >  0,  z  >  0 

(2) 

o 

II 

(3) 

was  originally  established  for  the  cases  where  the  matrix  M  €  /£"*n  is  either  a  P-matrix  (all 
principal  submatrices  have  positive  determinants)  or  else  is  PSD  (positive  semi-definite). 
See  [2], [3], [6], [5].  Furthermore,  in  all  these  papers  it  is  assumed  that  the  problem  at  hand 
is  nondegenerate,  meaning  that  for  each  solution  of  (1),  at  most  n  of  the  2n  variables 
u>i, . . . ,  wn,  2i, . . .  ,z„  are  zero.  Degeneracy  was  handled  by  allusion  to  perturbation  and 
lexicographic  techniques.  For  LCPs  with  P-matrices,  a  different  approach  was  proposed  by 
Murty  [10]  who  developed  a  Bard-type  algorithm  with  a  least-index  pivot  selection  rule.  The 
present  paper  extends  a  technical  report  by  Chang  [1]  in  which  degeneracy  problems  arising 
in  the  PPM  and  Lemke’s  method  [9]  are  resolved  by  means  of  least-index  pivot  selection 
criteria.  It  should  be  noted  that  the  results  of  [1]  are  confined  to  the  two  aforementioned 
matrix  classes. 

Recently,  Cottle,  Pang  and  Venkateswaran  [6]  introduced  the  class  of  sufficient  matrices  that 
contains  all  P-matrices,  all  PSD-matrices  and  other  matrices  as  well.  Actually,  the  class  of 
sufficient  matrices  is  the  intersection  of  two  other  classes:  the  row  sufficient  matrices  and  the 
column  sufficient  matrices.  As  shown  in  [6],  the  latter  two  classes  are  intrinsically  related  to 
the  linear  complementarity  problem.  Row  sufficient  matrices  are  linked  to  the  existence  of 
solutions  whereas  column  sufficient  matrices  are  associated  with  the  convexity  of  the  solution 
set.  (See  [6]  for  details.) 
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More  recently  still,  Cottle  [4]  showed  that  the  PPM  applies  to  nondegenerate  LCPs  with  row 
sufficient  matrices.  This  left  open  the  question  of  whether  a  least-index  degeneracy  resolution 
scheme  could  be  developed  for  this  extension  of  the  algorithm.  In  the  present  paper,  we  show 
that  when  the  matrix  M  of  the  problem  ( q ,  M)  is  (both  row  and  column)  sufficient  and  the 
least-index  pivot  selection  rule  introduced  by  Chang  (op.  cit.)  is  used,  the  PPM  will  process 
the  LCP  (9,  M)  in  a  finite  number  of  steps. 

Section  2  of  this  paper  contains  a  quick  review  of  some  background  needed  for  am  appreciation 
of  the  maun  result.  This  involves  the  basic  definitions  of  row  and  column  sufficient  matrices, 
a  review  of  principal  pivoting,  and  some  results  on  principal  pivoting  in  systems  with  row 
auid/or  column  sufficient  matrices.  In  Section  3,  we  state  The  Principal  Pivoting  Method 
with  least-index  pivot  selection  criterion,  and  in  Section  4  we  prove  that  it  processes  the 
LCP  (9,  M)  whenever  M  is  a  sufficient  matrix. 

2.  Background 

Since  the  matrix  classes  with  which  we  deal  are  not  well  known,  we  recadl  what  is  meant  by 
row  (and  column)  sufficient  matrices. 

Definition.  The  matrix  M  €  iT1*"  is 

(i)  row  sufficient  if 

Xi(MTx)i  <  0  for  all  t  =  1, . . . ,  n  =»•  x,(MTx),  =  0  for  all  *  =  1, . . . ,  n,  (4) 

(ii)  column  sufficient  if 

Xi(Mx)i  <  0  for  all  t  =  l,...,n  ==>  x,(Afx),-  =  0  for  all  *  =  1,. . . ,n,  (5) 

(iii)  sufficient  if  it  is  row  and  column  sufficient. 

Another  way  to  express  the  conditions  above  is  through  the  notion  of  the  Hadamard  product 
of  vectors  (or  matrices).  If  u  €  -ft”  and  v  6  iT*,  their  Hadamard  product  is  the  vector 
u*v  €  FC1  defined  by 

(u  *  v)i  =  Ui  • Vi  i  =  1, . . . ,  n. 

To  apply  this  notion  to  the  definition  of  a  column  sufficient  matrix,  we  let  u  =  x  and  v  =  Mx. 
Then  the  defining  condition  is 

x  *  (Mx)  <  0  =>  x  *  (Mx)  =  0. 

In  the  case  of  a  row  sufficient  matrix,  the  defining  condition  is 

x  *  (Mtx)  <  0  =>  x  *  (Mtx)  =  0. 

It  is  demonstrated  in  [4]  that  sufficient  matrices  are  different  from  P-matrices,  adequate 
matrices  (see  Ingleton  [8]),  and  PSD-matrices;  moreover,  they  are  not  necessarily  positively 
scaled  versions  of  such  matrices. 
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As  preparation  for  the  brief  summary  to  follow  and  for  the  statement  of  the  least-index 
principal  pivoting  method,  it  will  be  helpful  at  this  point  to  review  the  notion  of  principal 
pivoting.  In  equation  (1)  we  have  an  affine  transformation  of  n-space  into  itself  given  by 
z  H4  w  =  q+Mz  where  M  €  RnXn  and  q  €  K*.  For  the  present,  let  M  be  an  arbitrary  square 
matrix  with  the  property  that  for  some  index  set  a  C  {1, . . . ,  n}  the  principal  submatrix  Maa 
is  nonsingular.  We  may  assume  that  the  corresponding  principal  submatrix  of  M,  namely 
Maa ,  is  a  leading  principal  submatrix  of  M.  This  is  not  a  restrictive  assumption,  as  it  can  be 
brought  about  by  the  process  called  principal  rearrangement ,  the  simultaneous  permutation 
of  row  and  column  indices.  Now  consider  the  equation  w  =  q  +  Mz  in  partitioned  form: 

u>a  =  qa- f*  Maaza  -f  Ma&zs 

(6) 

ws  —  q&  +  MSaza  +  M&&Z8 

In  this  representation,  the  0- variables  are  nonbasic  (independent)  and  the  variables  are 
basic  (dependent). 

Since  Maa  is  nonsingular  by  hypothesis,  we  may  exchange  the  roles  of  wa  and  za  thereby 
obtaining  a  system  of  the  form 


za  =  q'Q  +  M'aawa  +  M'a&z& 
t"<s  =  <fs  +  MLwc  +  M^sza 

where 

=  -M£q.  AC  =  AC  AC  =  —AC  AC 

q&  —  Af<5Q Af00  qa  Mfio  —  Mg,aMaa  =  Mso ^,0 


Definition.  The  system  (7)  is  said  to  be  obtained  from  (6)  by  a  principal  pivotal  transfor¬ 
mation  on  the  matrix  Maa.  In  this  process,  the  matrix  Maa  is  called  the  pivot  block. 

The  following  facts  are  noteworthy. 

1.  Every  row  (column)  sufficient  matrix  has  nonnegative  principal  minors.  See  [6]. 

2.  Every  principal  rearrangement  of  a  row  (column)  sufficient  matrix  is  row  (column) 
sufficient.  See  [4]. 

3.  Every  principal  submatrix  of  a  row  (column)  sufficient  matrix  is  row  (column)  sufficient. 
See  [4]. 

4.  The  class  of  row  sufficient  matrices  and  the  class  of  column  sufficient  matrices  are 
invariant  under  principal  pivoting.  See  [4]. 


3.  The  Least-Index  Principal  Pivoting  Method 

In  this  section  we  tersely  state  the  (symmetric)  principal  pivoting  method  least-index  pivot 
selection  rule.  This  algorithm  (without  the  least-index  tie-breaking  rules)  has  previously  been 
extended  to  nondegenerate  linear  complementarity  problems  (q,  M)  in  which  the  matrix  M 
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is  row  sufficient.  The  proof  is  given  in  [4]  which  may  also  be  consulted  for  further  details  and 
references. 

The  PPM  works  with  pivotal  transforms  of  the  system 

tv  =  q  +  Mz.  (9) 

In  the  development  below,  we  use  the  superscript  v  as  an  iteration  counter.  The  initial  value 
of  v  will  be  0,  and  the  system  shown  in  (9)  will  be  written  as 

w°  =  q°  +  M°z°.  (10) 

In  general,  after  v  principal  pivots,  the  system  will  be 

wv  =  qv  +  Mvzv.  (11) 

Generically,  the  vectors  wu  and  zu,  which  represent  the  system’s  basic  and  nonbasic  variables, 
respectively,  may  each  be  composed  of  w  and  z  variables.  Principal  rearrangements  can  be 
used  to  make  {luj',  zj'}  =  {to,-,  z*}  i  =  1, . . .  ,n. 

Systems  like  (11)  are  traditionally  represented  by  a  tableau  (or  schema) 


The  symmetric  version  of  the  PPM  executes  principal  pivotal  transformations  (with  pivot 
blocks  of  order  1  or  2)  in  order  to  achieve  one  of  two  possible  terminal  sign  configurations  in 
the  tableau.  The  first  is  a  nonnegative  “constant  column”,  that  is,  gf  >  0  for  all  i  =  1, . . .  ,n. 
The  other  is  a  row  of  the  form 

q,  <  0  and  <0  j  =  1 , . . . ,  n. 

The  first  sign  configuration  signals  that  ( wv,zv )  =  (qy,0)  solves  (q,  M).  The  second  sign 
configuration  indicates  that  the  problem  has  no  feasible  solution.  The  PPM  (as  originally 
conceived)  does  not  actually  check  for  this  condition.  It  cannot  occur  when  mvTT  >  0  as  in 
the  case  of  a  P-matrix.  In  the  more  general  row  sufficient  case,  it  can  be  inferred  from  the 
condition 

q"  <  0,  mvrr  —  0  and  mvir  >0  Vi  ^  r, 
which  would  be  detected  during  the  “minimum  ratio  test.” 

The  PPM  consists  of  a  sequence  of  major  cycles,  each  of  which  begins  with  the  selection  of  a 
distinguished  variable  whose  value  is  currently  negative.  That  variable  remains  the  one  and 
only  distinguished  variable  throughout  the  major  cycle.  The  object  during  the  major  cycle 
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is  to  make  the  value  of  the  distinguished  variable  increase  to  zero,  if  possible.  Each  iteration 
involves  the  increase  of  a  nonbasic  variable  in  an  effort  to  drive  the  distinguished  variable 
up  to  zero.  This  increasing  nonbasic  variable  is  called  the  driving  variable.  According  to 
the  rules  of  the  method,  all  variables  whose  values  are  currently  nonnegative  must  remain 
so.  The  initial  trial  solution  is  (tw°,  z°)  =  (g°,0),  hence  at  least  n  of  the  variables  must  be 
nonnegative.  For  those  variables  whose  initial  value  is  qf  <  0,  we  impose1  a  negative  lower 
bound  A  where 


A  <  min  {of}. 

Then,  in  addition  to  requiring  all  variables  with  currently  nonnegative  values  to  remain  so, 
the  PPM  also  demands  that  the  variables  currently  having  a  negative  value  remain  at  least 
as  large  as  A.  This  broadens  the  notion  of  basic  solution ;  nonbasic  variables  are  now  allowed 
to  have  the  value  0  or  A.  (See  [1],  [2],  [3].) 

To  distinguish  between  the  names  of  variables  and  their  particular  values,  we  use  bars  over 
the  generic  variable  names  w"  and  zj'.  At  the  beginning  of  a  major  cycle  in  which  negative 
lower  bounds  A  are  in  use,  we  will  have  z?  =  0  or  z\  =  A  i  =  1, . . .  ,n.  Next,  we  use  the 
notation 

Wv{zv)  =  qv  +  Mvzu. 

The  definition  of  the  mapping  Wv  is  the  same  as  that  of  wv,  but  it  emphasizes  the  argument 
z". 

If,  at  the  outset  of  a  major  cycle,  the  selected  distinguished  variable  is  basic,  the  first  driving 
variable  is  the  complement  of  the  distinguished  variable.  Thus,  if  w VT  is  the  distinguished 
variable  for  the  current  major  cycle,  then  zvT  is  the  first  driving  variable.  The  distinguished 
variable  need  not  be  a  basic  variable,  however.  With  the  broader  definition  of  basic  solution 
(given  above),  the  current  solution  (u>*',z*')  may  have  zvT  =  A  <  0  at  the  beginning  of  a  major 
cycle.  In  such  circumstances,  zvr  can  be  the  distinguished  variable  as  well  as  the  driving 
variable.  In  this  event,  the  increase  of  the  driving  variable  will  always  be  blocked,  either 
when  a  basic  variable  blocks  the  driving  variable,  i.e.,  reaches  its  (current)  lower  bound  (0  or 
A)  or  when  the  distinguished  variable  increases  to  zero  (in  which  case  the  major  cycle  ends). 
The  point  of  the  least-index  degeneracy  resolution  scheme  presented  here  is  that  ties  in  the 
choice  of  the  blocking  variable  can  be  broken  so  as  to  insure  the  finiteness  of  the  PPM. 


1This  artifice  is  not  needed  when  it  is  known  that  M  €  P. 
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The  Least-Index  Rule.  In  applying  the  PPM  to  the  linear  complementarity  problem 
( q ,  A/),  break  ties  among  the  blocking  variables  as  follows: 

(A)  If  the  distinguished  variable  is  among  the  tied  blocking  variables,  choose  it  as  the 
blocking  variable  (and  terminate  the  major  cycle). 

(B)  Otherwise,  choose  the  (basic)  blocking  variable  with  the  smallest  index  as  the  exiting 
variable 

Symmetric  PPM  with  Least-Index  Pivot  Selection  Rule 

Step  0.  Set  v  =  0;  define  (tZ>°,  z°)  =  (q°,0).  Let  A  be  any  number  less  than  min; 

Step  1.  If  qv  >  0  or  if  («)*',  zv)  >  (0,0),  stop;  (tD*',  zv)  :=  (qu, 0)  is  a  solution.  Otherwise2, 
determine  an  index  r  such  that  zrr  =  A  or  (if  none  such  exist)  an  index  r  such  that 
wvT<0. 

Step  2.  Let  Cr  be  the  largest  value  of  zvT  >  z”  satisfying  the  following  conditions: 

(i)  zvT  <  0  if  zvT  =  A. 

(ii)  W?(z?,  •  •  • ,  *->,*?,  *n, -..,^)<0if<<  0. 

(iii)  W? (z . . . ,  zJLj,  z*+1,  •  •  •  ,  K)  >  0  ^  >  0. 

Ov)  wr (*? . . . . ,  «,  ^ i,  ...,«)>  a  a  t»r  <  o. 

Step  3.  If  Cr  =  +oo,  stop.  No  feasible  solution  exists.  If  Cr  =  0,  let  z"+1  =  0,  zf+1  =  Z*  for 
all  i  ^  r,  and  let 

^+i  =  v^‘'+1(r+1)  =  W'*'(z*'+1). 

Return  to  Step  1  with  v  replaced  by  v  +  1.  If  0  <  Cr  <  +oo,  let  a  be  the  unique 
index  determined  in  Step  2  by  the  conditions  (ii),  (iii),  (iv)  and  the  least-index  rule. 

Step  4.  If  m**1  >  0,  perform  the  principal  pivot  (xvvt,zvt).  Let 

z;+1  =  Wv,{z\, . . . ,  Zr-i,Cr,  Kn  •  •  •  ,*n)  and  w "+1  =  (r+1). 

If  s  =  r,  return  to  Step  1  with  v  replaced  v  +  1.  If  a  ^  r,  return  to  Step  2  with 
v  replaced  v  +  1.  If  mt'tt  =  0,  perform  the  principal  pivot  {(u£,z,),  (w* , zvt)}.  Put 
u>Jf+1  =  z*,t£*+1  =  Cr  5^+1  =  ^  f°r  *11  *  £  {r,s},  and  then  =  W-'+1(z‘'+l)  for 
all  i  ^  {r,s}.  Return  to  Step  2  with  z/  replaced  by  v  +  1  and  r  replaced  by  a. 


3 At  the  beginning  of  a  major  cycle,  for  each  index  r,  at  most  one  of  w* ,  z"  can  be  negative. 


4.  The  Finiteness  Argument 

In  the  following,  we  show  that  the  PPM  with  the  least-index  rule  (stated  above)  will  process 
any  linear  complementarity  problem  (9,  M)  in  which  M  is  sufficient.  It  is  interesting  to 
observe  that  the  mechanics  of  the  algorithm  itself  appears  to  require  only  the  row  sufficiency 
property.  The  finite  termination  of  the  algorithm  (with  or  without  the  least-index  rule)  is 
assured  if  each  major  cycle  is  finite,  for  the  total  number  of  negative  variables  is  nonincreasing 
during  each  major  cycle  and  decreases  strictly  at  the  end  of  the  major  cycle.  Such  finiteness 
is  realized  when  the  problem  is  nondegenerate.  We  show  here  that,  even  for  degenerate 
problems,  the  major  cycles  of  the  PPM  are  finite  provided  the  least-index  rule  is  enforced 
in  the  pivot  selection  criterion.  As  will  be  seen  below,  the  finiteness  of  the  PPM  with  the 
least-index  rule  hinges  on  the  column  sufficiency  property.  This  is  why  we  assume  the  matrix 
is  both  row  and  column  sufficient. 

If  cycling  occurs  in  a  major  cycle  it  is  not  restrictive  to  assume  that  it  is  one  in  which  u>i  is  the 
distinguished  variable.  It  follows  from  [4,  Theorem  4]  that  since  tt>x  and  z\  are  monotonically 
increasing,  both  u>i  and  Z\  are  fixed  during  cycling.  However,  the  algorithm  tries  to  increase 
Wx  or  Z\  in  this  major  cycle.  Hence  stalling  occurs  during  these  steps.  Accordingly,  if  we 
delete  all  the  variables  that  are  not  involved  during  cycling,  the  PPM  with  the  least-index 
rule  merely  looks  for  the  index  *  such  that 

s  =  min{*  :  <  0} 

and  then  pivots  on  muat  (if  mvat  0)  or  it  pivots  on 


Without  loss  of  generality,  we  may  assume  that  all  the  variables  are  involved  in  the  pivoting 
during  cycling.  Then,  during  cycling,  the  PPM  with  with  least-index  rule  performs  the  same 
pivoting  sequence  as  the  following  scheme  does. 

Scheme 

Step  0.  Start  with  the  system  tw"  =  qv  -(-  Mvzv ,  v  =  0,  where  w°  =  q°  +  M°z°  is  the 
initial  system.  (In  the  following,  M,”  represents  the  column  of  Mu  corresponding 
to  the  nonbasic  variable  z\  at  iteration  v.  Similarly,  M-'.  represents  the  row  of  Mv 
corresponding  to  the  basic  variable  u>{\) 

Step  1.  If  M.  1  >  0,  stop.  The  driving  variable  z\  can  be  increased  strictly.  Otherwise,  let 
s  =  min{t  :  m ft  <  0}. 

Step  2.  If  mvt,  >  0,  perform  a  pivot  on  mvtt  and  return  to  Step  1  with  v  replaced  by  v  +  1. 
Otherwise,  perform  a  block  pivot  of  order  2  on  the  principal  submatrix 

/  mu  mu  \ 

\  <1  ) 

and  return  to  Step  1  with  v  replaced  by  v  +  1 . 
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If  we  can  show  that  M.\  >  0  after  a  finite  number  of  pivots  in  the  above  scheme,  then,  since 
the  driving  variable  z\  can  be  increased  strictly  at  this  step,  we  obtain  a  contradiction  to  the 
assumption  that  cycling  occurs  in  a  major  cycle  (in  which  tvj  is  the  distinguished  variable) 
of  the  PPM  with  the  least-index  rule. 


Before  proceeding,  we  present  a  small  result  on  sufficient  matrices. 


Lemma  1.  Let  M  €  i?nxn  be  column  (row)  sufficient.  Then  for  any  real  numbers  a,  b,  c  such 
that  ab  <  0  <  c,  the  matrix 


/  rrin  mi2 
77121  77222 


M  = 


m„i  mn2 

\  b  o 


min  a  > 

m2„  0 

777 nn  0 

0  c  j 


is  also  column  (row)  sufficient. 

Proof.  It  suffices  to  prove  the  assertion  for  column  sufficient  matrices.  Let  x  =  (xx,  x2,  •  •  • ,  xn,  arn+i)T 
satisfy  the  inequalities  Xj(Mx)j  <  0  for  *  =  1, . . .  ,n  +  1.  Then  in  particular, 


xi(mnxi  +  . . .  +  rnmXn)  <  -axixn+i 


and 

bxiin+i  <  -cxl+i  <  0. 

Since  ab  <  0  it  follows  that  — axixn+i  <  0.  Thus 

n 

nijXj)  <0  t  =  1 , . . . ,  n. 
i= i 

Since  M  is  column  sufficient, 

n 

*  =  l,...,7i. 

;=i 

In  p«ui;icular,  it  follows  that  Xixn+i  =  0.  This,  in  turn,  implies  that  x  *  (Mx)  =  0.  □ 


Notice  that  Lemma  1  provides  a  mechanism  for  generating  sufficient  matrices  of  arbitrarily 
large  order. 

Lemma  2.  In  the  above  scheme,  a  pivot  in  row  s,  where  s  >  2,  must  be  followed  by  a  pivot 
in  some  row  with  a  larger  index  before  another  pivot  in  row  s  can  occur. 

Proof.  The  proof  is  by  induction.  If  the  matrix  M  is  of  order  1  or  2,  the  lemma  is  trivial. 
Suppose  the  lemma  holds  when  the  order  of  M  is  less  than  n  and  now  consider  the  case  when 
M  is  of  order  n. 
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We  shall  examine  the  situation  where  two  pivots  occur  in  row  a  and  2  <  s  <  n  —  1 .  If ,  between 
these  two  pivots,  there  is  no  pivot  in  some  row  with  a  larger  index,  then  by  deleting  M.n 
and  M„.,  a  contradiction  to  the  inductive  hypothesis  can  be  derived.  Therefore,  it  suffices  to 
show  that  there  is  at  most  one  pivot  in  row  n. 

Suppose  a  pivot  occurs  in  row  n  at  iteration  v\.  Let  (Tl)  denote  the  corresponding  tableau 
at  this  iteration.3 


1 

Z\ 

Ml 

9i 

mu 

*  •  •  rnln 

M„ 

9n 

mn  i 

^nn 

Tableau  (Tl) 


By  the  choice  of  the  pivot  row,  we  have  mfl  >  0  for  all  i  <  n  —  1  and  m„x  <  0  in  (Tl). 

Suppose  the  next  occurrence  of  a  pivot  in  row  n  is  at  iteration  i/2.  When  this  occurs,  zn  must 
be  the  exiting  basic  variable  and  Wi  is  either  basic  (Case  I)  or  nonbasic  (Case  II). 

Case  I.  (tox  is  a  basic  variable  at  iteration  v2.)  Let  <r  be  the  set  of  indices  i  such  that  w,  is 
nonbasic  at  iteration  i/2.  Note  that  1  £  a.  Let  M  denote  the  principal  transform  of  M  at  this 
iteration.  Clearly,  M  can  be  obtained  from  M  by  performing  a  block  pivot  on  the  principal 
submatrix  Maa.  Thus, 


Max  =  -M£M9X. 

Now 


(12) 


©  N 

f  ©  \ 

0 

© 

. 

and  Ma\  — 

I 

© 

© 

L  —  1 

L  —  1 

(13) 


Being  a  (nonsingular)  principal  submatrix  of  a  sufficient  matrix,  M„„  is  also  a  sufficient 
matrix.  (See  [4].)  From  (12)  and  (13),  we  have 


f  ©  > 

( e\ 

(  ©\ 

© 

e 

© 

*  (ACM,,)  ~ 

: 

* 

: 

! 

© 

© 

© 

\  “  ) 

<  +  > 

l  -  y 

3For  simplicity,  we  represent  (Tl)  without  using  superscripts. 
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which  is  impossible  since  Ma)  is  column  sufficient. 

Case  II.  (u^  is  a  nonbasic  variable  at  iteration  i/2.)  Let  the  definition  of  er  be  as  in  Case  I,  but 
note  that  now  we  have  1  €  a.  Since  Maa  is  sufficient,  the  diagonal  entry  mn  is  nonnegative. 
There  are  two  cases. 


Case  II.  1  (mn  >  0.)  The  pivot  on  mn  would  not  change  the  sign  configuration  of  MaX 
namely 


©  \ 
0 

© 

\  “ 


Once  this  pivot  is  performed,  we  have  Case  I  (with  a  different  index  set  a). 


Case  n.2  (mn  =  0.)  Here  there  are  two  more  cases. 

Case  II.2.1  (mn  >  0.)  By  performing  a  pivot  on  mn  in  schema  (Tl),  the  variable  wx 
becomes  nonbasic  and  the  sign  configuration  of  M.x  is  unchanged.  Therefore,  as  in  Case  I,  a 
contradiction  can  be  derived. 


Case  n.2.2  (mn  =  0.)  Let  (T2)  denote  the  tableau  at  iteration  v2. 

*<T 

Tableau  (T2) 


1 _ w„  z9 


9^ 

A/<7<7 

Mg# 

9* 

Mg(f 

M99 

In  this  tableau,  a  =  {1, . . .  ,n}  \  a. 

Now  let  qn+x  €  R  be  arbitrary  and  enlarge  (Tl)  to  (Tl*)  as  follows 


1 

Z\ 

z2 

zn 

zn+l 

wx 

9l 

mn 

mX2 

•  •  •  min 

-1 

U)2 

92 

m2i 

m22 

•  •  •  m2n 

0 

U>n 

9n 

mnX 

mn2 

•  •  ’  mnn 

0 

«>n+l 

9n+l 

1 

0 

0 

1 

Tableau  (Tl*) 
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By  Lemma  1,  the  bordered  matrix  of  tableau  (Tl*)  is  sufficient.  The  block  pivot  on  the 
principal  submatrix  Maa  in  (Tl*)  produces  a  tableau  (T2*)  having  (T2)  as  a  subtableau. 


1 

W<j 

to 

Zn+1 

to 

to 

Maa 

Afj.n+l 

to 

Mgg 

Mgtn+1 

tfn+1 

9n+l 

Mn+ltg 

■^fn+l,n+l 

Tableau  (T2*) 


Notice  that  (T2*)  has  the  same  basic  ^-variables  as  (T2),  hence  Tableau  (T2*)  is  the  corre¬ 
sponding  enlargement  of  (T2). 


Pivot  on  Moa 

(Tl) 

(T2) 

Enlargement 

i 

|  Enlargement 

(Tl*) 

(T2*) 

Pivot  on  Moa 

By  pivotal  algebra,  we  have 

^hn+1,1  ~  (Afn+i^Mffo 

)i 

—  Mn+ 1  ,<T  i 

=  (l,0,...,0)Mffl 


=  mn 

=  0. 

Now  the  matrix 

mll 

^n+1,1  ^n+l.n+l 

is  nonsingular  and  can  be  used  as  a  pivot  block  in  (Tl*).  Denote  the  resulting  tableau  by 
(T2**).  In  it,  tui  and  tun+1  are  nonbasic  while  all  the  other  u>,  are  basic. 
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1 

U>1 

*2 

zn 

«>n+l 

Z\ 

mn 

mu 

•  •  •  mi„ 

mi(n+i 

W2 

h 

m2i 

rhn 

•  •  •  m2„ 

m2,„+i 

Wn 

9« 

m„i 

mn2 

mnn 

m„,„+i 

Zn+1 

9n+l 

m„+i.i 

m„+i.2 

mn+i,n 

mn+i,n+i 

Tableau  (T2**) 


In  Tableau  (Tl*),  we  have 

mu  =  0,  mu  >0  i  =  2,...,n  —  1,  m„i  <  0,  and  m„+iti  =  1. 


Hence 


mn  =  1,  m,i  >  0  t  =  2, . . .  ,n  —  1,  m„i  <  0,  and  m„+i,i  =  -1. 

Now,  since  both  (T2*)  and  (T2**)  are  principal  transforms  of  tableau  (Tl*),  it  follows  that 
(T2**)  is  a  principal  transform  of  (T2*).  In  fact,  if  we  define  the  index  set  p  —  (cr  \  {1})  U 
{n  + 1},  then  (T2**)  can  be  obtained  by  performing  a  block  pivot  on  the  principal  submatrix 
M„p  in  (T2*).  Therefore  Mpi  =  —M~^MpX.  The  indices  n  and  n  +  1  belong  to  p  and 


m„i  <  0,  mn+i,i  =  0,  m„i  <  0,  mn+,t,  =  -1 
while  m,i  >  0  and  m,i  >  0  for  all  other  i  €  p-  Accordingly,  we  obtain 


f  ©  > 

i  e\ 

/  Q\ 

Afpi  *  (M~pMpi)  ~ 

0 

* 

© 

© 

— 

+ 

— 

1  o  y 

(  +  y 

^  0  y 

But  this  is  impossible  since  M~p  is  (column)  sufficient.  □ 

Lemma  3.  In  Tableau  (Tl),  M,\  >  0  after  a  finite  number  of  iterations. 

Proof.  For  j  >  1,  let  p(j)  be  the  number  of  pivots  that  occur  in  row  j.  In  the  proof  of 
Lemma  2,  we  have  shown  that  p(n)  <  1.  Furthermore,  it  follows  from  Lemma  2  that 

M  <  Z)  M  0  +  1- 

«'=>+! 


12 


In  other  words, 


fi(n  —  1)  <  fi(n )  +  1  <  2 
n(n  —  2)  <  22 

fi(n  -i)<  2*-1  +  2‘-2  +  •  •  ■  +  2  +  2°  =  2*. 

Therefore,  the  scheme  will  terminate  after  a  finite  number  of  iterations.  □ 

Theorem.  In  the  case  of  a  linear  complementarity  problem  (q,  M)  with  a  sufficient  matrix 
M,  every  major  cycle  of  the  PPM  with  the  least-index  rule  consists  of  a  finite  number  of 
pivots. 

Proof.  Suppose  cycling  occurs  in  a  major  cycle  in  which  wx  is  the  distinguished  variable. 
Then,  since  wx  and  zx  are  monotonically  increasing,  both  wx  and  zx  are  fixed  during  cycling. 
However,  it  follows  from  Lemma  3  that  M.x  >  0  after  a  finite  number  of  steps.  Therefore 
either  wx  or  z\  can  be  strictly  increased  after  a  finite  number  of  steps,  in  contradiction  to  the 
assumption  that  cycling  occurs.  □ 

Corollary.  In  the  sufficient  matrix  case,  the  PPM  with  least-index  rule  will  process  the  LCP 
( q ,  M )  in  a  finite  number  of  steps. 

Proof.  Each  major  cycle  of  the  algorithm  reduces  the  number  of  negative  components  in 
(w,  z)  by  at  least  one.  The  assertion  now  follows  from  the  Theorem.  □ 

Remark.  In  implementing  the  least-index  rule  it  is  important  to  obey  statement  (A)  which 
says  that  if  the  distinguished  variable  is  among  the  tied  blocking  variables,  then  it  is  to  be 
chosen  as  the  blocking  variable.  Failure  to  do  so  can  lead  to  the  false  impression  that  the 
problem  is  infeasible. 
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